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Abstract 



We analyse the Langevin dynamics of the random walk, the scalar field, 
the X-Y model and the spinoidal decomposition. We study the deviations 
from the equilibrium dynamics theorems (FDT and homogeneity), the asymp- 
totic behaviour of the systems and the aging phenomena. We compare the 
results with the dynamical behaviour of (random) spin-glass mean-field mod- 
els. 



1 Introduction 



Spin-glasses and other disordered systems have 'critical' dynamics through- 
out their low-temperature phase. Their most striking dynamical effect is that 
of aging: they do not reach thermal equilibrium after very long times and 
experiments are performed out of equilibrium showing a dependence on the 
history of the system |^ . 

These phenomena have been studied with numerical simulations P, ^ ^ 
P] , several phenomenological models have been proposed [0, H , and analytical 
techniques have applied to mean-field systems |P, |10|, [ill . 

It has been recently proposed that mean-field spin-glasses do not to 
reach a situation of dynamical equilibrium {i.e. homogeinity in time and the 
fluctuation-dissipation theorem (FDT) are violated) even after infinitely long 
times. 

In other systems an intermediate situation known as interrupted aging 
occurs in which non-equilibrium effects tend to dissapear, but very slowly as 
compared to the relaxation of ordinary non-critical systems {e.g. paramag- 



nets, ferromagnets) ||T2|, [13 



The persistence of out-of equilibrium effects after very long times, and 
in particular the violation of the equilibrium theorems is a feature not re- 
stricted to disordered systems such as spin-glasses. It is interesting in itself to 
study the deviations from the equilibrium theorems in simpler examples with 
hamiltonians that are deterministic (non random), and even not disordered 
or frustrated. 

The aim of this note is to analyse the Langevin dynamics of some such 
simple examples: the random walk, the free scalar field and the X-Y model. 
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and the spinoidal decomposition of a ferromagnetic Ising-like system. In 
each case we shall study the deviations from the equilibrium theorems and 
we shall analyse the long-time behaviour of the correlation and response 
functions, and the total response to a constant perturbation applied during 
a finite time-interval (the equivalent of the 'thermoremanent magnetization' 
in spin-glass experiments). 

The organisation of the paper is the following. In section 2 we present 
some general remarks on the FDT and its possible generalization. In section 
3 we analize the simplest non equilibrium model, i.e. the random walk. In 
section 4 we consider the case of the D-dimensional free scalar field theory. In 
section 5 we study the dynamics of the X-Y model at low temperature in two 
dimensions. Finally in section 6 we consider the dynamics of the spinoidal 
decomposition for the usual ferromagnetic Ising case. Our conclusions are 
presented in section 7. 



2 The generalized fluctuation dissipation re- 
lation 



Let us consider a system which has been quenched from high temperature at 
time t = 0. The auto-correlation function C{t,t') among a local quantity O 
at two subsequent times t' and t is 




(2.1) 



Hereafter ( • ) represents the mean over the thermal noise. 
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For large t and fixed t — t', in an equilibrium dynamics situation the 
auto-correlation function behaves as 



C{t,t')=C{t-t') , (2.2) 

i.e. it is homogeneous in time. 

The response function to a pertubation is defined as the variation of the 
quantity {0{t)) with respect to a perturbation applied at time t'. More 
precisely, if we consider the perturbed Hamiltonian 

H = Ho + J dt h{t)0{t) , (2.3) 

the response function is defined as 

and, because of causality, it equals zero if t' > t. The response function is not 
independent of the correlation if the system is in equilibrium. Indeed, it is 
related to the correlation function by the celebrated fluctuation-dissipation 
theorem (FDT): 

R{t,t') = M{t-t')^-^^, (2.5) 

and it is also homogeneous in time -R(t, t') = R{t — t'). 

However, if the system is out of equilibrium neither homogeneity nor the 
FDT ( |2.5| ) hold. The generalised relation between response and correlation 
functions can be written as 

RM=m-t')X{t,t')^-^^^ (2.6) 

with X a function of both times t' and t that characterizes the approach to 
equilibrium. 
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The FDT and its violation can be partially understood from the following 
considerations. Let us consider a system described a variable y{t) which 
satisfies the Langevin equation 

j^y{t) = -F[y]{t)+r^{t) (2.7) 

where 77 is a Gaussian random noise with zero mean and correlation 

{v{t)r^{t')) = 2T6{t-t'), (2.8) 

T being the temperature. 

Taking t > t' for definiteness, the equation of motion ( |2.7D implies 

(I7 - |) C{t, t') = 2TR{t, t') + A{t, t') (2.9) 

where we have used ( y{t) rjit') ) = 2T R(t, t') and 

A{t,t') ^ {F[y]{t)y{t')-F[y]{t')y{t)) . (2.10) 

At equilibrium the correlation functions satisfy ( B(t)D{t') ) = { B{t')D{t) ), 
if B{t) and D{t') are any two functions of y(t). This is a consequence 
of the time reversal symmetry. Hence the asymmetry A vanishes and the 
fluctuation-dissipation theorem may be recovered by using the invariance 
under translations in time of the correlation functions at equilibrium: 

C{t,t') = C{t-t') (|, + |)c(t,t') = (2.11) 

and 

R{t,t')=f3^^. (2.12) 



In the off-equilibrium situation the homogeneity in time (eq.( p.ll| )) is not 



valid and the asymmetry A may be present. Eq. (2^) is not valid in general 
and the generalisation (|2.6|) must be considered. 



In mean-field spin-glass models the auto-correlation and response func- 
tions are defined as C(t,t') = (1/A^) Si(t)si(t') ) and 
R{t,t') = (l/A^)E^i'^(s»(t))/5/ii(t'), respectively. In the analysis of the 
asymptotic dynamics presented in Ref. [0 (see also [Q) it has been pro- 
posed that, for long enough times and small time differences, t,t' ^ oo and 
{t — t')/t << 1, X = 1 and FDT is satisfied, while for long enough times and 
big time differences, t,t' ^ oo and {t — t')/t ~ 0(1), the function X depends 
on the times only through the correlation function C(t, t'), i.e. 

R{t, t') = (39{t - t') X[C{t, t')] . (2.13) 

A self-consistent asymptotic solution for the mean-field out of equilibrium 
dynamics within this assumption has been found both for the p-spin spherical 
and the Sherrington-Kirkpatrick models (Refs. |9|1,|[T0 . 



In the following sections we shall investigate the behaviour of the function 
X for various (non random) models and we shall compare the results with 
expression ( |2.13| ) at long times. 



2.1 Scalings and aging phenomena 



Another interesting problem is to study the scaling properties of the correla- 
tion, response and X functions, and the response of the system to a constant 
perturbation applied during a finite period [0,t^]. 

^Asymptotic means t^t' ^ oo after N oo 
This result is expected to hold for large times for systems with finite susceptibility. 
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For the p-spin spherical spin-glass model, if t and t — t' are both large one 
analytically finds 

(2.1.1) 

within the assumptions described above. The numerical solution of the dy- 
namical equations suggests that 

C{t,t') = S{t'/t) , (2.1.2) 

[i.e. ha. power law). This is a new (non- homogeneous) scaling Q This is the 
simplest scaling that captures an essential feature of spin-glass phenomenol- 
ogy: the aging effects, i.e. the explicit dependence of the behaviour of the 
system on its history. The scaling ( p.l.2|) can be modified in many ways 
to describe in more detail the results of simulations of realistic models and 
experiments. Then, one sometimes assumes the slightly different form 

C{t,t') =t-^ S{t'/t) , (2.1.3) 

where 5 is a small number, of the order of a few percent in spin-glass models 
[0, P, The factor implies an interruption of aging for large t (0(few 
years)). 

The generalised FDT relation ( |2.6| ) can be written as 



R{t,t')=pQ{t-t')X[t-\t'/t] -^C{t,t'), (2.1.4) 



If we substitute the scaling form ( |2.1.3| ) for the correlation function in eq. ( |2.9| ) 



and we assume that a similar form is valid for the asymmetry, we find that 

the response function scales as C/t' or dC/dt' . 

^In the dynamics of other mean-field spin-glass models more complicated scalings can 
be present. 
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If instead we assume that the asymmetry A is zero, as will be the case 
below, we find 

Xit,t')=XiX) = ^ 

with A = t'/t. 

In the typical aging experiments 0] one measures the 'thermoremanent 
magnetization', i.e. the response of the system to a constant magnetic field h 
applied during the interval [0, t^,], at constant temperature, is interpreted 
as a 'waiting time'. In a general dynamical system described by the Langevin 
equation ( |2.7D the equivalent of the thermoremanent magnetization is 

XtAt)= / dt'Rit.t'). (2.1.6) 

Jo 

Aging experiments show that Xt^ if) = (t) depends non-trivially on t and 

tn. Si- 



l + X + 6 



S'{X) 



(2.1.5) 



3 The Random Walk 



The simplest example of a dynamical system that does not reach equilibrium 
is the random walk |T^. In the continuum limit the quantity y{t) satisfies 
the very simple differential equation 

j^yf) = r/(t), (3.1) 

with T] a Gaussian noise with variance given by eq.(^]8|). 
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It is easy to check that the correlation function C{t, t') = { y(t) y(t') ) and 
the response function R{t,t') = S{yit) )/6h(t') = {13/2) {y{t)rj{t') ) are given 
by 

C{t,t') = 2T min(t,t') , (3.2) 
R{t,t') = 6{t-t'). (3.3) 

Hence, the relation (|2.6|) is satisfied with 

X(t,t') = l/2, (3.4) 

Vt,t', a constant function but different from the usual FDT result, X = 1, 
the system never reaches equilibrium. 

Indeed one finds that the scaling form (|2.1.3 ) for the correlation is satisfied 



with S{X) = A for A < 1, but with a big value for 6, 6 = —1. Inserting 
this scaling in eg. ( p.1.5] ) we also obtain X = 1/2, as expected since in the 



random walk problem the force F and the asymmetry A are zero. However, 
the scalings for the correlation and the total response are quite different 
from those observed in spin-glasses. In terms of the 'waiting time' ty^ and 
T = t -t^, C(r + t^,t^) = 2Ttw and Xt^i^ + tw) = tyj {cfr. eq. (|2.1.6|) ). 
Both expressions are independent of r but depend explicitly on tw 

This example may seem trivial, but it captures the essence of the phe- 
nomenon that will be described in the rest of this note. 
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4 Free Gaussian Fields 



In this section we study the behaviour of a simple free scalar field 0(x;t). 
The Hamiltonian is quadratic in the field and in dimensions D it reads 



{W<PY + mi<P' (4.1) 



where mo is the mass of the field (see e.g. Ref . [ITTj ) 



The relaxational dynamics is given by the Langevin equation 
d 

— 0(xo; to) = A0(xo; Q - ml (f){xo] Q + r]{xo; to) . (4.2) 

?7(xo; to) is a Gaussian noise (^^(ko; to) its Fourier transform) with zero mean 
and correlations 

{r^i^o;to)vK;t'o)) = 2T exp(-2^) 5(t„-0 
{viK; to) viK'; to)) = 2T i2n)^ exp{-§) 5^ (K + K) 5 (to- t'o) 

Xo = |xo — x'qP and = |kop. We have introduced a spatial correlation 
over a typical lenght 1/A to simulate the lattice spacing. (This serves to cure 
some large k pathologies.) 

Taking 0(xo, 0) = as the initial condition, the solution to the dynamical 
equation ( [4. 2] ) for each noise realisation is 

0(xo; to) = J ^ e^^-° dr e-(^^+™?)(*-^) r^{K, r) . (4.3) 

Since we are dealing with a field, the correlation, response and X functions 
depend on space-time coordinates. A standard calculation for the correlation 
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function Co(xo, x^; to, t'J = { 0(xo; to) t'^) ) gives 

The response function -R(xo, Xq'; ^o) — ^( 0(^o; ^o) ) /9h{^g] t'^) is given 

by 



(47r)^/2 {to -t'^ + ^ri^ • ^ 

Here and in what follows we take unprimed times bigger than primed times 
and we omit the theta functions. 

The preceeding formulae suggest to measure space, time and mass in 
appropriate lattice units 

t = AHo, (4.6) 

X = Ax„, (4.7) 

m = rrio/A (4.8) 

and to rescale the correlation and response functions C = A^'^~^^Co, and 
R = ( 1/A^ )i?o. Note that there is no rescaling of fields and correlations in 
D = 2. 

The function X that measures the departure from FDT reads, in terms 
of the rescaled coordinates: 
Xfx x'- 1 t') = y -R(x,x^;t,0 

dt' 



(t-t' + l\'"^ [ . 2/ ^'i' 

1 + exp <^ -2mH' + — 

\t + t'+l) ^\ 2[(t+l)2-t'2] 



(4.9) 
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4.1 Large-times behaviour 



Consider first the massive case. We have a time scale given by: 

teg ~ . (4.1.1) 

For any x fixed and any two times t,t' » t^q we have that X = 1. This 
identifies teq as an 'equihbration' time |. 

In the massless case m = the equihbration time diverges and we have a 
more interesting situation. Let us concentrate in this case. For fixed x and 
large times t, t', we have 



X{x]t,t') = X{x;\) = 




(4.1.2) 



with A = t'/ 1. Hence X is non-trivial and FDT is violated, even for very long 
times. 

If A ^ then 

X(x;A)^l/2. (4.1.3) 

If A ^ 1 and L> ^ 0, then 

X(x;A)^l, (4.1.4) 

and we recover FDT. A = 1 corresponds to times t, t' satisfying {t—t')/t « 1, 
i.e. small time differences. 

If we put = we recover X = 1/2 for all times, the result for the 
random walk. 



^Note however that if x is of order ^/iJ or larger then X can be smaher than 1, even 
zero for small time differences. We shall not consider such diverging distances in the rest 
of the section. 
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4.2 Scalings 



We now present the scalings. Since the massless scalar field turned out to 
be more interesting we shall concentrate in this model. If m = the explicit 
computation of the integrals in eq. (|4.4|) gives 



C{^,^';t,t') = T , r 



X 



X 



D 

~2~ '4(t + t' + l)'4(t-f + 1) 



where r[?T,; a, h] is the generalised incomplete Gamma function 



r [n; a, h] 



(4.2.5) 



(4.2.6) 



For equal space points x = eq. (|4.2.5|) reduces to 
1 1 



C(0;t,t') =T 



(47r)^/2 1 



D 



(t + t' + 1)1-^/2 + 1)1-^/2 _ 

(4.2.7) 

For long times and A < 1, i.e. big time differences (t — t')/t ~ 0(1), this 
expression satisfies the scaling law ( 2.1.3| ) with 5 = D /2 — 1 and 

1 



S{\) = T 



^^.(|l_A|i-^/2 + |i + A|i-^/2) 



(4.2.8) 



(4vr)^/2 (1 - f ) 

Hence, in this time scale the function X (eq. ( [4.1.2| )) can be written as 
X(0; t. A) = X[t^C] and in particular X(0; t. A) = X(C) for D = 2. 



Considering again the general model, the total response (|2.1.6|) reads 



c/^x / dt" 



-^(t-*")-Z(i477TT) 



(477)^/2 {t - t" + 1)^/2 



-m?{t-tw) 



1-e 



-m''tii 



(4.2.9) 
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and for large t^, >> teq it reduces to 



1 



— m r 



(4.2.10) 



the typical relaxation in a system that has equilibrated; i.e. Xt^i'T') depends 
only on T = t — tw and no aging is present. 
Instead, in the massless limit 



which shows a dependence on the history for all t^, although a rather unusual 
one. 

The learned reader will notice that in the massless case the Hamiltonian 
is invariant under the transformation 



The correlation functions are not invariant under this transformation and 
therefore the symmetry is spontaneously broken. The slow approach to equi- 
librium is a reflection in the time domain of the Goldstone boson arising from 
the spontaneous breaking of the symmetry. In the next section we shall see 
a case where the symmetry group is implemented in a non linear way. 

These results are in agreement with the general formulae discussed in 
Section 2 when the asymmetry is neglected (c/r. eq. ( |2.1.5| )). Indeed the 
asymmetry is zero, because the force F is linear in the field 0. 




(4.2.11) 




(4.2.12) 
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5 The relaxational dynamics of the XY model 



The Hamiltonian of the 0(2) non-hnear a model can be written in terms of 
the angular variable defined through S(xo) = (cos 6'(xo), sin ^(xq)). In two 
dimensions it reads 

(5.1) 



H = -j d'^o {ve{^o,Qy 

(see e.g. Ref. ||T7[). 

The relaxational dynamics is given by the Langevin equation 

— ^(Xo; Q = - + r/(Xo; Q , 

oto ot^(Xo; to) 



(5.2) 



with //(xo; to) as in eq. (|4.3| ). 

We consider low temperatures such that vortices can be neglected and 
therefore we do not see the Kosterlitz-Thouless transition. 

The solution to the dynamical equation (|5.2| ) for each noise realisation is 
that of the massless scalar field problem in D = 2, cfr. eq. (|4.3|) . 

The angle-angle correlation C{^Xo, x^; to, t'^) = { 6{x'^; t'^) 6'(x„; to) ) is given 

by 



T 

c{xo] to, t^) = — r 



2^2 



0; 



2^2 



(5.3) 



4(l + A2(t„ + t^))'4(l + A2(t„-t^))_ 

and in particular, the time correlation between the angles at the same space 
point (Xg = 0) is 



An ^\l + A\to-t'J 



(5.4) 



The response to an external field h{'Xo;to) acting like h{'Xo;to) 9{^o'-,to), 

i?„(x„,x;;t„o = 6{e{^o;to))/6h{^',;t'^) = i/{2T){e{^o;to) v{<;t'o)) is 
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given by 

R(r-f f')- ''''^^ ^(l+AHto-t'o))) (rr) 

Ko[Xo,to,tJ - l+A2(t,-t^) ■ ^^■^> 

We now turn to calculating the physical quantities for which the angular 
character of 9 is essential. We first calculate the 'composite' correlation 

C{xo;to,0 = (sin^(x„;g sin^(x;;0) (5.6) 
= exp (^-i (^(0; to, to) + ^(0; t'„ t',))^ sinh C{xo; to, t',) 

(5.7) 

and the associated response to a transverse field /i(xo; to) acting like 
h{y.o,to) sin6'(xo;to) 

where m(xo;io) is the transverse magnetisation m{^o',to) — ( sin ^(xo; to) )■ 
The 'composite' response can also be written in terms of the angle-angle 
correlation C and its associated response function Ro'. 

Ro{xo;to,t'J = e^p(^~{C{0;to,to) + C{OX,t',))+C{xo;to,t',)^ 

X Roixo;to,t'^) . (5.9) 

As in the previous section we now rescale space-time coordinates as 

t = AHo , (5.10) 
X = Axo , (5.11) 

and rescale the response function R= {l/K^)Ro, R= {l/h?)Ro (but neither 
the correlations nor the angles). 

In terms of the new coordinates we have 



C{x;t,t') = £r 



0; 



4(l + t + t')'4(l + t-t') 



(5.12) 
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R{x;t,t') 



1 



1 



An 1+t-t 



7 exp 



A{l + t-t') J ' 

C{x;t,t') = {{l + 2t){l + 2t')y^ smhC{x;t,t') , 
R{x;t,t') = ^{{l + 2t){l + 2t')y^ R{x;t,t')exp(^C{x;t,t')) .(5.15) 



From eqs. ( ^.14| ) and ( |5.15| ), the function X reads 

2R{x;t,t') 



X{x;t,t') 



A^{x; t, t') + A+{x; t, t') exp{-2C{x; t, t')) 



with 

A±{x;t,t') 



dC{x;t,t') ^ l dC{0;t',t') 
dt' 2 dt' 



T 

47r 



exp 



A{l+t-t')) ^ ( 4(l+t+t')) _^ 1 



(5.13) 
(5.14) 



(5.16) 



l + t-f 



l + t + t' 



l + 2t' 



(5.17) 



5.1 Large-times behaviour 



In this subsection we consider the large times hmit, t and t' large [t > t'). In 
this limit the function X is 

lim Xix- 1, t') = lim ^^-^ ^ (5.1.1) 

i^oo *^°°l + exp(^-2C{x;t,t')j 

with C{x; t, t') given by eq. (|5.12|) . We shall analyse the function X and the 
correlation C in different regions determined by the space and time separa- 
tions X and t — t' . 
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Equal Times 

We first consider the correlation and response functions at equal times t = 
t' ^ 1. We consider separately the cases x = (local values) and x ^ 1 
(many 'lattice spacings'). 

a. X = 

C{0;t,t)^^ (5.1.2) 

which was to be expected, since the 0{2) symmetry is unbroken and 
( sin^ 9) = 1/2. We also have 

X{0;t,t)^l (5.1.3) 
i.e. the system evolves locally with an equilibrium dynamics. 

b. > 1 

b.i t>a;2>l 

lim X{x;t,t) = 1 

t—*oo 

and 

C{x]t,t) Cstatic{x) ^ X~^ 



(5.1.4) 
(5.1.5) 



b.ii a;^ > t > 1 

\im X{x;t,t) = J, 

t— >oo Z 

C{x;t,t) 0. (5.1.6) 

We conclude that if we take a snapshot of the system at a large time t, 
within a range of length ~ t^/^ the system seems equilibrated in the sense that 
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the correlations coincide with the static ones and the response satisfies 



FDT. Well outside that range the angles are uncorrelated and X = 1/2, as 
in a random walk. In the following we shall see in more detail the nature of 
this 'equilibration'. 

Different Times 

We here consider different times, i.e. t — t' ^ and we again analyse sepa- 
rately the cases x = and x ^ 1. 



a. X = 



lim X{0;t,t') = (5.1.7) 

' ' t+t' \"2;r 



i+t-t' 



and the correlation C{0;t,t') reads 



T 
47r 



CiO-,t,t')^li,uri{{j^y-(^] 1 (5.1.8) 



b. Consider now two points widely separated x ^ 1 but well within the 
'range of equilibration' for these times, i.e. 

t > t' > > 1 (5.1.9) 

Two possibilities then arise: 
b.l {t-t')/{t + t') 

(A — > 1.) We reobtain the 'equilibrium' situation X = 1 and the correlation 
function goes, asymptotically in t — t', as 

C{x;t,f) ^ {t-t')~^ . (5.1.10) 
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b.tt {t -t')/{t + t') > 

(A<1.) 



lim X{x]t,t') = X{\) 



1 



(5.1.11) 



and the correlation function is 



C{x;t,t') 



(2t)-£ 
2 



A »^ 




1 - av- 
i + aJ 



(5.1.12) 



(c/r. eq. (^X8|) ), and X = t-^/(4")X[C]. 



For relatively small time differences with respect to the total time and for 
any two points well within a 'domain' of equilibration these results are still 
those of a system evolving as in equilibrium, in other words X = 1 and the 
correlation and response functions are homogeneous in time. However, even 
within a 'domain' when the time separation is large enough, the correlation 
and response give manifestly out of equilibrium results: X has a non-trivial 
time dependence and the correlation and response functions are not necessar- 
ily homogeneous in time. We conclude that one cannot picture these domains 
as regions in which a true (lasting) equilibrium has been established. This 
has to be contrasted with the behaviour of the massive scalar field, which 
after a certain teq and for fixed x evolves as in equilibrium. 
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5.2 Remanent magnetization 



Let us now turn to studying the behaviour of the system in a 'thermorema- 
nent magnetization' experiment. 

The total response to a constant, uniform magnetic field h{x; t, t') applied 
from if:' = to t' = tyj over the whole system is 

xtjt) = l(i + 2t)"^ fj dt' [1 + 21')-^ ^(tTTtI'^) ^^-^-^^ 

with 

^) = ^ du exp (^-|^ exp (-^r 
Defining a function (see Appendix A) 

/,..r)^.-.^....-.4^(l^)-^Kr,, (5.2.3) 

we have that in the large t limit for every and t 

XtJt)=tl-^f{j,T) , (5.2.4) 

and f{Xw, T) is an increasing function of X^, finite at = 1. Asymptotically, 
for t':^tuj 

Xt^t) (X tl'^ (5.2.5) 

Several considerations are in order about this behaviour. Throughout the 
low temperature phase T < Air, and the susceptibility diverges with tw. This 
result was to be expected since the static magnetization grows as h'' with 
7 = T/(87r-T) 0. 



u uw 



(5.2.2) 
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For finite times and small fields the linear response theory holds, but 
becomes worse as an approximation for larger times and it fails completely 
at oo. This result is reminiscent of what seems to happen in spin-glasses 
with the reaction of the system when temperature is slightly changed: while 
experimentally (long times) this response is possibly non-symmetrical with 
respect to the sign of the temperature changes, it is still symmetrical in the 
relatively short times involved in most simulations. 

The longer the waiting time during which the field has been applied, the 
slower the relaxation of the remanent magnetization. If we normalize the 
magnetization by its value at ty^, then the decay is a function of t/tw This 
is again reminiscent of what happens in spin glasses and other disordered 
systems, except for the fact that the susceptibility is finite in those cases. 



6 Spinoidal Decomposition 



We consider a normal ferromagnetic system (of Ising or Heisenberg type) 
and we supose that the dynamics is local, without local conservation of the 
magnetization. The Langevin equation described in the previous sections is 
a good example of such a dynamics. For definiteness we consider the Ising 
case. 

We are interested in studying the evolution of such a system when we 
quench it from high temperature to a subcritical temperature. The problem 



is well studied in the literature |T^. The main result is the random formation 
of domains oriented in different directions which become larger and larger 
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with increasing time. The size of the domains ^(t) grows as t^^^. It is also 
well known that the equal time correlation function in the large time limit is 
well described by 

(<^(x,t)(/.(0,t))=F(x/e(t)), (6.1) 

for well separated space points; i.e. x » 1. The function F is not very far 
from a Gaussian. 

Throughout this section the brackets stand for average over initial condi- 
tions. We suppose that the field at time zero is Gaussian-distributed with 
a correlation function that goes to zero at large distances. 

The most natural proposal for the correlation function at different times 

is 

{(f>{x,t')(j){0,t)) = C{x/at),t/t') (6.2) 

and, if we consider the correlation at the same space point we would then 
have 

{mt')mt)) = c{t'/t) . (6.3) 

Intuitively we can understand this scaling as follows. The correlation 
function is proportional to the probability that both spins stay in the same 
cluster. At time t the spin stays in a cluster of size ^{t) which has a mean 
life proportional to ^(t)^, i.e. to t. Therefore it takes a time of order t to 
revert the magnetization. 

Let us be more precise. We consider the following zero temperature 
Langevin equation: 

^0 = A0 + 0(l-02). (6.4) 



One can treat this problem approximately as follows [|T9|: introduce a 
field ip defined by 
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Eq. ( |6.4| ) becomes, in terms of ip, 

The approximation consists in neglecting the last term; then one assumes 

^ = A^ + ,p. (6.7) 

Note the sign of the mass term. We do not discuss here the range of validity 
of this approximation, which is widely done in the literature flSl . 



The strategy we follow is similar to the one used in the preceding sec- 
tion: we first solve a simple linear problem and then calculate the physical 
correlations as correlations of composite operators. 

The ^/'-correlation function in Fourier space is given by 

(^(k,t)^(k',t')) oc5^(k + k') exp{{t + t'){l-k^)) (6.8) 

which corresponds to 

( ^(x, t) ^(x', t') ) oc (^^^ exp(t + t' - ^^^) (6.9) 

in position space. For large times the absolute value of iJj becomes exponen- 
tially large, goes to ±1 (c/r. eq. (|6.5|)) and one has 

C(t, t') = ( sgn(7A(0, t)) sgn(V>(0, t')) ) • (6-10) 

The correlation of the random Gaussian variables ipi = ifj{:x. = 0,t) and 
i/j2 = ipi^ = 0,t'), is 

i^^^j)^ ^t^^\)D/2 ^''^'' ^,J = 1,2, (6.11) 
their randomness comes from that of the initial conditions. 
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Eq. ( |6.1CI| ) becomes 



J dij, #2 exp [-j^- - 2(t> j 

Changing variables 

exp(tO 

- f^^^ (6-13) 

we obtain 



r(f f'] = sgn(^i)sgn(^2) exp {-Ajjl - 2Bipi^2 - ^2) . 

^ ' ^ / #1 #2 exp (-A V^f - 2B ^1^2 - V^i) ' ^ ^ ^ 

with 

A = X-I^ B={^-^y\ (6.15) 

A = t' /t. Since we are only interested in the scaling we do not explicitly 
compute this integral, but notice that 

C{tX)=C{X). (6.16) 

In a similar way one can prove the other scaling laws. The important 
result is that this approximation gives expressions for the correlations which 
are in very nice agreement with the aging formula . 

A similar analysis could have been done for the response function and 
for the correlation at finite temperature, but this would make this paper too 
long. 
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7 Conclusions 



We have seen that in many systems in which equihbrium is slowly approached 
some form of aging phenomena are present. These systems are characterized 
by a correlation length that is infinite in the static limit but is finite for finite 
times: it diverges with a power law in time. 

A remarkable feature is that in these systems the energy landscape is 
fiat: no high barriers in energy are present. On the contrary the flatness of 
the potential in certain directions, i.e. the presence of zero modes, is at the 
origin of this very slow approach to equilibrium. 

These systems are an evident proof that it is not possible to conclude 
for the existence of energy activated barrier-crossing only from the presence 
of aging. It would be rather interesting to see if there are some peculiar 
phenomena, which may distinguish the effects of barriers from those due to 
fiat directions. 



Appendix A 



The total response to a constant magnetic field /i(x; t) applied during the 
interval [0,^^] over the whole system defined in eq. ( p.l.6|) is 

XtAt) = \{l + 2t)~^ J^^dt'{l + 2tr^ k(^\±^,T^ (A.l) 

with 

k{w,T) = j du exp (-^ ) exp (-— r 



u uw 



(A.2) 
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Changing variables 



(A.3) 



the total response is 



l + 2(l+t) 



1 -CJ 

1 + u; 



k{u,T) 
(A.4) 



Let us now consider the A scale, i.e. t ^ 1 and < At„ = t^/t < 1. In 
this case we can use 1 + at ~ t and 1 + 1 — t' ~ t — t' in the lower limit of the 
integral. Then, 



du 



l + \v 



1 — cuX 



and finally 



k{uj,T) (A.5) 



(A.6) 



with /(A^,t) given by eq.( |5.2.3| ). Note that the integrand in ( |A.5|) diverges 
for u; = 1 and gor uj = 0, but the integral over u is still convergent, and ( |A.5| ) 



is valid for all A^. 
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